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ated manuscript. 
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AN ANALYSIS OF THE SHAPES OF A CELL DURING DIVISION 
WITH PARTICULAR REFERENCE TO THE 
ROLE OF SURFACE TENSION 


H. D. LANDAHL 
THE UNIVERSITY OF CHICAGO 


An expression is given for the contribution to the relative elongation 
by any volume element or its surface. The case in which only a constant 
surface tension is acting is considered in greater detail. From measure- 
ments of an egg in several stages, the empirical variation of an expres- 
sion proportional to the above function is obtained. In the constricting 
region and near the ends of the cells, the change of shape is opposite 
to that which would be due to a constant surface tension alone. The 
effect of streaming which may arise from a variable surface tension is 
considered. The effect of forces which arise in this manner may be 
sufficient to explain the discrepancy in the constricting region if the 
streaming is considered due largely to such a variable surface tension. 


The shape of a cell at any instant, as well as changes of shape 
with time, can be considered as the resultant of various volume and 
surface forces. If the former can be derived from a potential, one 
may substitute for them an equivalent surface force. The average 
relative change in length of the cell can be written as (Rashevsky, 


1940) 
ied NAB [eZ —4(2X + yY)] ds (1) 
é=— = 24.—4(% yY ; 
r, at av | 3 


where 2r, is the length, 7 is the viscosity, V the volume of the cell, 
X , Y and Z are the components of the total surface force given by the 
sum of the surface forces and the surface force equivalent of the vol- 
ume forces. If we consider a figure of revolution about the z-axis, let 
p? = x? + y’ and let v be the normal to the surface at (p, 2), then equa- 
tion (1) may be written, on integrating with respect to the polar 
angle, 

Sdemahe ae 255) sec (yp), (2) 

6z  37V 


where P is the component of the total surface force along the p-axis. 
The left hand member of (2) represents the contribution to aver- 
age relative elongation @, by the surface of the element between z 


and z + 6z. Its value can be measured for any z ata particular time, 
so that Z(P) could be determined. This gives little information about 
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Z and P,, which themselves are composite surface force components. 
We shall ipetend consider the component of 6é/é6z is due to a surface 
tension, y, which has a constant value over the surface. If R, and 
R, are the principle radii of curvature and R twice their harmonic 
average, then that part of dé/dz due to this factor may be written 


dé 2ny p 2ay 

— = — tan (vp) — Se &); 3 

is 3,0 R [z (vp) — 4p] 3yV > ) (3) 
where ¢(z) is defined by this equation and is a function of z since the 
shape determines p(z). The function ¢(z) can be determined empir- 
ically for any given shape by measuring p, R, , R. and tan (vp) for any 
z. The empirical expression for ¢(z) may give some insight into the 
role played by surface tension. For this reason, enlarged photographs 
of three successive stages of an Arbacia egg, denuded by treatment 
with KCl, were measured to determine $(z) as shown in Figure 1. The 
photographs were obtained through the courtesy of Dr. R. Buchsbaum. 
Each point is obtained from measurements of the egg. The curves 
correspond to the tracing of the shapes correspondingly numbered in 
the figure. The time between stages is half a minute. The curves, 
which are free hand drawings through the points, are stopped at the 
end point of the cell along the elongation axis. The value of ¢ at the 
end of the cell is obtained from measurements and the limit of p tan 
(yp) is obtained analytically by assuming that the end is approximat- 
ed by part of a sphere. The area under the positive part of each curve 
is very nearly equal to the total area under the negative parts. From 
the curves, we see that the effect of a true surface tension for these 
shapes is such that in the region of the neck the tendency is to reduce 
elongation, and thus to reduce constriction. The effect becomes pro- 
gressively greater with time. Evidently some other surface or volume 
force than a surface tension must be acting here to cause constriction. 
In the middle range, from about .5<10-° to 3X10- em., the effect of 
surface tension is such as to produce an elongation. This effect, as 
measured by the area under the positive part of the curves, is much 
the same in successive shapes, perhaps decreasing with time. If so, 
the change in shape in this region is as expected from a surface ten- 
sion. The part of the curve near the end shows the existence of a ten- 
dency to reduce elongation. A comparison of the areas under the 
curves indicates that the tendency to oppose elongation becomes great- 
er with succeeding shapes. Thus the shape in this region is changing 
in the wrong way from that expected by a surface tension. Thus it is 
probable that the effect of volume forces and surface forces, other than 
that due to a constant surface tension, must be such that, in the neck 
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region and near the ends, the contribution towards elongation will be 
positive. 

The variation of surface tension from pole to equator has been 
suggested as a factor in cell division (Spek, 1918, cf. also Wilson, 
1925). It should be noted that a variation of the surface tension with 
2 would not greatly affect the above argument as such variations of 
force could be sustained by the outside medium. However, this is per- 
haps not probable. It is more likely that variations of surface tension 
would be sustained largely by streamings set up in the cell. In gen- 
eral, such streamings are observed. One might attempt to compute 
the magnitude of the velocity of the streaming necessary to sustain 
the variations required to make the net elongation positive, or, using 
the observed order of magnitude of such streamings, to obtain the 
effective volume force in the neck region and see whether this could 
overbalance the negative effect of surface tension. 


2 
z IN 10°CM. 
FIGURE 1 


Consider the case of the shape II of Figure 1. The average of ¢ 


in the neck region in which ¢ < 0 is ¢ = —1.4X10° cm. If we use 
+ © 10> dynes cm.+, the average value of 6é/éz, in this region, due to 
surface tension is —.5/n sec... The contribution to 6é/6z for this re- 
gion due to streaming, which we assume to be due only to the variable 
surface tension, may be estimated in the following manner. If the 
velocity of streaming in this region, is assumed to be ~ 10° cm. sec.”, 
falls off to a negligible value in a distance, d ~ 3X10~* cm., or about 
one-half the width of the neck, so that d= z for ¢ = 0, then the force 
per unit area is 7v/d. The area upon which this is acting is of the or- 
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der of zr2, where r is the radius of the neck. This force is equivalent, 
in order of magnitude, to a pressure on the surface of the neck, in the 
region where ¢ < 0, of ar? n v/(dard). Dividing this pressure by 37 
we obtain for the contribution to 6é/dz in this region due to the 
streaming, the value rv/(3d?)% 10 sec.? We have considered the 
pressure to be of the order of magnitude, of the pressure difference. 
Comparing this value with that for surface tension alone, and giving 
yn a value of 10+ poise, we see that the latter could be an important 
factor at this particular stage. In this illustration we may have un- 
derestimated the former and overestimated the latter. It may be noted 
that the effect of the viscosity does not enter the resulting expression 
for the contribution to the relative elongation rate due to streaming. 
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EQUILIBRIUM SHAPES IN NON-UNIFORM FIELDS 
OF CONCENTRATION 


H. D. LANDAHL 
THE UNIVERSITY OF CHICAGO 


The special case of a metabolizing cell, with infinite permeability, 
suspended in the field of an equivalent cell is treated exactly to deter- 
mine its equilibrium shape for small deformations. It is found that flat- 
tening takes place along the line of centers. 


In the early cleavages of an egg, the cleavage plane generally de- 
pends upon the relative position of the particular blastomere. On the 
basis of the diffusion force theory of cell division (Rashevsky, 1940), 
one would in general expect such a result. For the present we confine 
ourselves to a very idealized case of a cell of radius 7), which metabo- 
lizes a substance of molecular weight M at a rate q gm cm® sec", 
whose permeability to the substance is infinite and which contains 
structures. on which diffusion forces may act. We consider first 
the effect of the presence of a second cell upon the first cell when the 
second is similar to the first but not deformable, i.e., we neglect 
the effect of the first cell on the second. We take the center of the first 
cell as the origin, the line of centers as the z-axis and the distance 
between centers as d. The cell will be deformed at a rate which is 
given by the expression (Young, 1939) 


o= | [Cz cos(v 2) — 4[x cos(y xz) + ycos(y y)]} dS (1) 


where 
© = — RT ue/2n.MV;; (2) 


c is the concentration of the metabolite, 7, is the viscosity of the cell, 
V, its volume, and u is a measure of the relative volume of particles 
in the cell. Actually the expression for ® is more complicated than 
given above (Landahl, 1942a), and the above is correct only as to or- 
der of magnitude when the particles, upon which the force is acting, 
are small. 

In addition to the deforming force, there will be a repulsion be- 
tween the cells. We shall not consider this effect for the present. 
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Introducing spherical coordinates, 7, @ and @, and setting é = 
cos 6, equation (1) becomes 


e= | [ orsp.(edédy, Sg 
0 -1 


where P(é) is the second Legendre polynomial and 7’ is the value of 
y at the surface of the cell. We shall consider only very small deforma- 
tions, so that we may set 7” = 7%. 

The concentration due to the first cell at P can be expressed by 


seid : 


where D, is the external diffusion coefficient, p, is the distance from 
the element dV to the point P , and the index 1 refers to the first cell. 
Since we are considering only very small deformations from the 
spherical shape, we have 


6, (P) =¢1 Vi/4aD ar « (5) 
Similarly for the second cell we have 
63(P) = q2V3/4nD 15 (6) 


where 7, is the distance from the center of the second cell to the point 
P. Butif H(é,7/d) is the generating function for the Legendre poly- 
nomials, then (6) may be written 
Tor 
Q2V2H (é ’ dae 


C2 (P) ad Wide adaet (7) 


Introducing ¢ = ¢, + c, into (2) and (8), we obtain 


-E "fle tHe, a | Oeste, (8) 


Yor 
or 
2) eat 277 01 No Ke 
;  eiBeny 2 pe 
where Rashevsky (1940) 
K=RT qur*,/4MnD.. (10) 


If y, is the surface tension of cell one, we have the contribution to 
€ , due to this force, the amount (Young, 1939) 
fe V1 (11 — Ts) 


rs ae i @1) 
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where a coefficient which varies from slightly greater than unity to 
slightly less, as 7,/72. varies from 1 to o, is neglected. Thus setting 
the sum of these terms equal to zero we have the equilibrium defor- 
ae Here 27, and 2r, are the length and width of deformed cell 


e= (1% — To) /Tor (12) 
we have, on neglecting powers of « greater than the first, 
Ee — 8 Ne Ton Ke/15 V1 a3 . (13) 


Thus since ¢ is negative, cell one is compressed in the direction of the 
line of centers. 

If the parameters of the cells are identical, the maximum value 
of €, €m, is given when d = 27, or 


Em = — Ky/1dy. (14) 


For dividing Arbacia eggs Ky/y turns out to be about 2/3 (Landahl, 
1942b). Thus «,, is of the order of five per cent for this case. For eggs 
not about to divide the value of Ky/y has not been estimated, but if 
this turned out to be about one per cent or less, it would be very diffi- 
cult to detect such an effect because of difficulties of measurement. 
However, a much smaller effect might be ample to determine the di- 
rection in which the cell will elongate if this effect takes place at a 
proper stage. Thus generally the cleavage of two daughter cells is in 
the direction as required by the above considerations. 

If we pass to the case in which there are four cells whose centers 
form a square, each cell touching two others, we can determine that 
the plane of cleavage would be in the plane of the cells, as is generally 
the case. 

From the method of images we would predict that there should be 
a flattening effect, in addition to gravity, when a cell rests upon an im- 
permeable plane, while a cell at the intersection of two planes should 
elongate in the direction of the line of intersection of the planes. Such 
effects could show up in the predominance of the orientation of the 
cleavage planes in the proper directions. 

In the above development, we have assumed an infinite permea- 
bility, and thus the effect from one cell to the next may be greatly 
different in magnitude from that predicted above, though in the same 
direction. For this reason, the estimates of magnitude are little more 
than illustrations. 

In view of equation (3), we may determine é for rather general 
fields of concentration. Let the field be independent of y. The concen- 
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tration, and thus © of equation (2), is given by a sum of spherical 
harmonics. If @ is thus expressed, and A, and B, are coefficients of the 
harmonics of degree n and (—n-—1), then we have from equation (3), 
e=2V,(A: to, + Be 7701) (15) 
for very small deformations. At equilibrium the deformation is then 
obtained in a manner similar to that used in deriving equation (18). 
Thus, for small deformations, a knowledge of A. and B, is sufficient to 
determine the effect of the concentration field. 
This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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CELLULAR FORMS: THE TRI-AXIAL CELL. I 


ALSTON S. HOUSEHOLDER 
THE UNIVERSITY OF CHICAGO 


Without assuming a spheroidal shape for the cell, more general ap- 
proximate deformation equations are derived which allow for non-uni- 
formities in the external diffusion field and within the cell itself. It is 
found that in the absence of such non-uniformities, however, the only 
stable shapes are those in which the cell is spheroidal. Moreover, in con- 
formity with previous results derived by somewhat different approxima- 
tions, it is found that the spherical equilibrium is always stable and the 
approach monotonic unless, for a sufficiently large volume or a sufficient- 
ly high metabolic rate, there exist two non-spherical equilibria, one of 
prolate type and one of oblate type. 


1. The diffusion equation. Rashevsky’s discussion of diffusion 
in the cell, of cell division, and, more generally, of cellular forms 
(1940, chapters i, iii and v) presupposes the cell to be roughly sphe- 
roidal in shape. He thus ascertains the conditions under which the 
cell will be either oblate, prolate, or spherical. The present object is, 
essentially, to carry out the same calculations without, however, mak- 
ing the spheroidal assumption, and so to determine the more general 
equilibrium shape, i.e., its diameter along each of three perpendicular 
directions. The further possibility is also allowed that the external 
diffusion field has been rendered unsymmetric—say by neighboring 
cells unsymmetrically placed, or by other means. This is a possibility 
considered by Williamson in a forthcoming paper, but with Rashev- 
sky’s assumptions of symmetry. 

In attempting to approximate in the usual fashion we must revise 
somewhat the details of the geometric picture. Instead of approximat- 
ing the cell shape by a spheroid or by a cylinder with spherical caps, 
we must now approximate it by a rectangular parallelepiped of di- 
mensions 27, X 27, X 2r,. There are six faces to the cell and we may 
denote by the subscript pairs (71) and (72) (j =1, 2, 3) the two 
faces normal to the 7-th axis. We may then, as usual, let cj, and Cje 
be the concentrations just inside, c’;, and c’;. the concentrations just 
outside the corresponding: faces. 

In deriving the approximate diffusion equations Rashevsky in- 
troduces the quantity 6 to characterize the external diffusion field. 
This is the distance from the cell membrane at which, if the concen- 
tration gradient were constant, the concentration would fall to the 
value c, which it has at infinity. Rashevsky supposes 6 to be the same 
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on all sides of the cell. We shall remove this restriction. Moreover, 
Rashevsky supposes that the “average” value of the internal concen- 
tration is at a point halfway from center to periphery. We shall leave 
unspecified this point of average concentration and say only that 
aj7; OY 4,7; represents the distance from the center to the face (71) 
or (j2)-at which this average concentration occurs. With this under- 
standing we may write in the usual way 

C — Cio C sa, — Cy 


= h (jo ao Cia) =D. 


Ojo 1} af ’ (o LAR hs 2, 3), ( ) 


D; 


thus assuming the principal diffusion directions to be outward. 


If we set 
B=E = Op > (2) 


then we can solve the equations (1) for cjo and ¢jc: 


, + 1 + De nsen! Sy vj sje 
(ee (6 en ——— 
7 ( h djo0 D; Ojo : 


ae (1 ay ee )( a9 is Bi aja 7; De Ne 
Cio = —— : 
‘ : h Ojo h Ojo D; Ojo ; 


In each second there passes through the face (jc) the amount 


(3) 


iy D. a;¢7; De\ 
Anh (Cj — C jc) = 47 =(1 + sitet :) 
LUC; jo) K av ee a ae C5 (4) 
whereas the total production within the cell in one second amounts to 


8g 721s , and the increase in the amount of metabolite within the cell 
is 87,727; dc/dt. Hence 


(5) 


h dj0 h ajo YT; -1 
1 
: + os) re 


dc WSS FI 

— — ihe 
q 2 3( ts D: D; 

The coefficient of ¢ in this equation is the value of A-', the reciprocal 
of Rashevsky’s “diffusion resistance” for this case. For constant 7, 
the equation is again linear in ¢ and requires no discussion. We turn 
therefore to the deformation equation. 

2. The deformation equations. The derivation of these equa- 
tions (Rashevsky 1940, chap. iii) requires the evaluation of three sur- 
face integrals, each of six terms having the same form except for 
constant multipliers. The first integral, describing the effect of the 
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volume diffusion forces, consists of terms of the type 
Cia T; Sjo— 4r, 121; Cjo— 3 VEjo, 


with a constant multiplier. Since the entire integral involves also the 
terms 
— 1 V (Cro + Cec), 


it is legitimate, and convenient, to replace each ¢js by Cjo — C. 
The second integral describes the effect of the osmotic forces and 
a typical term is 


1; Xj0 Sjo, 


where Xjo is the component of the osmotic pressure. But this pres- 
sure has the direction of the normal and it is proportional, at the face 
is, 60 Cig — Cia — (Cj, — C) — (Cice—_Cg). The multiplier. V enters, 
just as before. Apart, then, from the multiplier (37V)+ which affects 
all the integrals, the combined first terms of the first two integrals 
amount to 


3RTu kT 
eaz V (Cia — Co). + oar (Cit aCe), aC 35.— Co) | 
RT 3 h dj0 hdjo, hoajo 1; \ 
om pessa(ee elie? D, aa D, ye, 


and the complete expression is obtained by summing as to o, and 
subtracting half the sum of the corresponding expressions obtained 
when the subscript 7 is replaced by k and by /. 

The third integral expresses the effect of the surface tension. 
Here a difficulty arises, for if we regard the cell as a parallelepiped 
the pressure due to the surface tension will be zero along the sides 
and infinite along the edges. Actually the curvature will be approxi-. 
mately that of an ellipsoid and we take 


ya aT) 
as the approximate pressure on either face S;, so that the force is 
3 7 Amy (Fyn tele V4 shes | 
and the corresponding term of the integral for the two opposite faces 
Ko —8y7; (41 + 7%) =— V y(t tr). 
This entire integral is therefore 


Vy (ri? — ore — 4"). 
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To combine these results, we have now found for the approximate 
deformation equations the following form: 


dr; RT 3 h dj0 h dj 
ry > nn) 1+ 


ot 
dt 67M 2 D. D. 
(6) 
h ajo 7; [7 ise ; 
ties cacti, snerekh Be 


where the round brackets {.--} denote the sum of the two expressions 
which are obtainable from the one written by changing the subscript 
j to k and to lin turn. These equations are derived from Betti’s for- 
mula for plastic flow, and this, in turn, presupposes constant volume. 
We note, indeed, that the sum of the right members of the three equa- 
tions (6) vanishes identically: 


Sy pt dr; =0 ’ 
which is equivalent to the statement 
1, To Ts; — const. 


3. Equilibrium states. If we denote by 6;? the quantity exhibit- 
ed on the right in (6), then these three equations have the form 


v5 

arg Oke a (Oe? =20)-9 (7) 
These quantities 6; have the dimension of time. Equilibrium occurs 
when the right members of (7) vanish, and this occurs only when the 
6’s are all equal: 


tee Pipe 


6, = 02 =6,;= 8, (8) 


say. In fact, if the 6’s and the a’s are regarded as constant (e.g. the 
a’s equal to 4 and the 6’s to half the cube root of the volume), we see 
that 


6; =6(153 Qj1 , jo » Oj: , Ojo). (9) 


That is, regarded as functions of the variable r and the four para- 
meters, the functions 6 have the same functional form. In principle, 
equilibrium may be determined as follows. Solve 


6(7; 01, 02,6,, 62) =6 (10) 
for 

r= R(0;0,,02,6;, 62), (11) 
Set a, = aj: , +++, do = dj to get 
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7; — K; (6), (12) 
solve 


8h, (6) R2(0) R3(6) =V (13) 
for 6, and finally insert this value of @ into (12). 


In this procedure the concentration ¢ enters as one of the para- 
meters and we have thus obtained the v’s as functions of c. If the 
variation of ¢ is regarded as negligible this completes the determina- 
tion of the equilibria. However, with variable ¢, we have obtained 
only the equilibrium values of the 7’s as functions of c. To complete 
the determination one must substitute these functions of ¢ into the 
right member of (5), equate this to zero, solve for c , and employ this 
value for the numerical calculation of the equilibrium 7’s. 

Unfortunately we cannot expect, in general, to obtain the equi- 
libria in closed form since the equations (9) are quadratic in the 7’s, 
the functions R thus involve radicals, and when (13) is rationalized 
we may expect an algebraic equation of rather high degree. We are 
forced, therefore, to consider solutions and approximations in special 
cases. 


4. Types of equilibria in a uniform field. When the external dif- 
fusion field is uniform and the cell is non-polar, we may at least drop 
the second subscripts on the 6’s and the a’s. Rashevsky, in fact, iden- 
tifies all the 6’s, and also all the a’s, and we shall do likewise. If, 
more generally, the differences of the 6’s , and those of the a’s, are to 
be regarded as vanishing only with the differences of the 7’s, second 
and higher approximations can be obtained from the present con- 
siderations by introducing and solving for small deviations. 

Rashevsky regards the 6’s as equal to the smallest r. It simplifies 
the discussion considerably if we regard the 6’s as equal to half the 
cube root of the volume of the cell, and this we shall do. Then drop- 
ping subscripts, the form of 6 is given by 


Aig eal, (14) 
where 
har iene vha Se 
oye Cite hp] Do) © paca ls (ene Dy eae” 
(15) 
RTD; 


3 
2—- —___|- 4 (1 + hp/D. -1] ; 
ape! oe 2 
The quantity represented by /? could be positive or negative, and for 
the moment we may disregard the positiveness implied by the expo- 
nent. 
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If we let 
hap 


7— 6A juic == a 
D,(1 + hp/De) 


(16) 


then the equilibria, for a given c, are given by the solutions 2; and 7 
of the four equations 


ris ast — #1 + 4;)>, 


(17) 
XH 3 —o". 
An obvious solution is the equilibrium of spherical type 
V4 = Ve — U3 — =>, (18) 


and we shall later consider the stability of this equilibrium. For the 
present we wish to inquire as to the existence and character of equi- 
libria of other types. 

Non-spherical equilibria are given by solutions of (17) in which 
at least two of the x’s are different. If x, # x. , then 


C0 = BLS a) a ee) 

or, after suppressing the factor x, — %2: 

Bate (ay Sad ea (19) 
Similarly, if +; # x2, 

Posti—~@. cH) (ee): 
But these two equations imply that 

ti = 23, 

whence no equilibrium is possible for a non-polar cell in a uniform 
field with the three dimensions all unequal. 


For the non-spherical equilibria, therefore, we have to solve equa- 
tions (19) and 


Ts 6 (20) 


for 2, and x. Now the graph of (19) is a rectangular hyperbola 
which intercepts the axes at (—1, 0) and at (0, —1), and for which 
the line «, = a, is always an axis of symmetry. Since (19) may also 
be written 
1 1 jie 
(23 — p) (te RD =(B=b?’ (21) 


we see that the center lies at [(6? — 1), (6? — 1)~]. If f? < 0 the 
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hyperbola does not intersect this fixed axis of symmetry at all, so that 
neither branch of the hyperbola enters the first quadrant and no in- 
tersections with (20) occur. Thus the equilibrium of spherical type is 
the only one possible. We therefore consider only the case of 6? > 0. 

If 6’ < 1 the hyperbola intersects the fixed axis of symmetry at 
XL, = X= (+8 — 1) in the third quadrant. Hence again the hyper- 
bola does not enter the first quadrant and only the spherical equilib- 
rium is possible. 

If 6? > 1 one branch is entirely in the first quadrant. It is geo- 
metrically evident that these curves have two intersections (possibly 
coincident) or else none, and that for a given o there are always two 
intersections if 6 is sufficiently large. If 


o > (B—1)2 (22) 


the curve (20) intersects the axis of symmetry of (19) inside the hy- 
perbola and hence the two non-spherical equilibria exist and one is of 
prolate type and. one of oblate type. If (22) fasls, there are either no 
intersections, or else there are two, both corresponding to equilibria 
of the prolate type. Two equilibria of the oblate type cannot exist. 
This conclusion is at variance with a result previously obtained 
(Householder, 1941) where, by extending Rashevsky’s calculations 
based upon the assumption of a spheroidal cell, conditions were found 
under which two equilibria of oblate type could exist. There 6 was 
taken as equal to 7, , which presumably accounts for the discrepancy. 

The relation (22) failing, there is a minimal value of §, for a 
given oc, for which the curves (19) and (20) intersect. This can be 
obtained by eliminating x, between the two equations and equating to 
zero the discriminant of the resulting equation. This yields a relation 
between f and co. Since both quantities involve the volume, the ex- 
pressions for f and o in terms of the volume can be substituted in this 
and we obtain thus an equation giving the minimal volume of a cell 
which is capable of having non-spherical equilibria, the internal con- 
centration being supposed fixed. The result of this calculation is 
omitted. 

Relation (22) may be written in the form 


pa (ere lee". (23) 


and if, in this, we insert for 6 and o their expressions in terms of p 
we obtain, when fractions are cleared and all terms transposed to the 
left, a cubic inequality in p. In this the leading coefficient is positive 
and the last two negative. Thus there is a minimal volume for which 
there exist non-spherical equilibria of opposite types. 
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We have so far neglected the variation of c. When this is con- 
sidered, the equilibrium is determined by equating to zero the right 
member of (5), an equation which may be written 


Q=P Dar + a), (24) 


h? 2? My 
Q = 0/.B eB 3 . (25) 
RTD;?(1 + hp/De)*15 a + hp/ De) 4h 


We need consider only the case 6? > 0, and when this holds we have 
B>O. Equation (24) together with (19), (20), and x, = #; give the 
required non-spherical equilibria, with 6”, which is proportional to c, 
regarded as one of the unknowns. We eliminate 6? from (24) by 
means of (19) and obtain 


ted 2(ay hl) 


aioe Tp Ws 


which, in view of (20), is equivalent to 


This is an ellipse whose eccentricity is 1/.\/2, whose major axis is 
horizontal, and which intercepts the x-axis at the origin and at 
v%,—=-—1. For Q (and hence q) sufficiently large the ellipse enters the 
first quadrant, and given any p a q may be determined sufficiently 
large so that (26) and (20) have two real intersections in this quad- 
rant. The intersections exist and correspond to equilibria of opposite 
type (cone oblate and one prolate) whenever 


Qe asl) (27) 
In case (27) fails, there is a narrow range of values of Q for which 
two non-spherical equilibria of the prolate type exist, while for all 


values of @ below this range the only equilibrium is the spherical one. 
5. Stability of the spherical equilibrium. Let 


x; =o0(1+ &;). (28) 
Equations (7) have then the form 
dé 
A(1+ gains os bi) ich pte hy) shea a 5) 
b? 
a Gales lar pall ge awit a — fy) (29) 


= CL giceiagah 
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If we expand in power series and observe that 


O—6,76 +6 4+---, (30) 
we have 
dé; Bo? (1a)? 3 
dt = aetige eo! oy 
where terms of degree higher than the first have been omitted. For 
a fixed value of c, therefore, the equilibrium is stable provided 


fo <A ss), (32) 


that is, provided either 6? is negative, or else 6? is positive and (22) 
fails. Otherwise the spherical equilibrium is unstable. Hence when 
non-spherical equilibria of opposite types (one prolate and one oblate) 
exist, the spherical equilibrium becomes unstable, but in all other cases 
the spherical equilibrium is stable and the approach to this equilib- 
rium is monotonic. 

To include the variation in the concentration, we have to adjoin 
equation (5), which can be written 


de 
Gq Pie Peed set p) ae ye (33) 
Since /? is proportional to c: 
ice 
say, the equilibrium value of c is 
«Vell to) - 
pee REY Een (34) 
and if we set 
C= 6, (1. &), (35) 
and expand, we have 
dé 3AB : (36) 
ane Seay es 


If we compare (25) and (15) and recall that A represents the co- 
efficient of c in (15), we note that the product AB is always positive. 
The four linearized equations to be considered are now the three 
equations (31) in which f? is evaluated with c = c, and equation 
(36). Evidently the stability of the equilibrium is governed only by 
the signs in (31), since that in (36) is always negative, so that the 
conclusions with regard to stability remain as given. 
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THE LINEAR THEORY OF NEURON NETWORKS: 
THE STATIC PROBLEM 


WALTER PITTS 
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_ The construction of a theory of activity in neuron networks of ar- 
bitrary topological structure is commenced under the linear excitation 
hypothesis: we consider conditions for possible steady-state equilibria, 
deferring a dynamical treatment to the sequel. 


I. Introduction and Preliminary Definitions. 

In four previous papers published in this Bulletin (1941a, b, ¢, 
1942; to be hereinafter referred to as I, II, III, IV respectively) A. S. 
Householder has considered the nature of steady-state activity in 
nerve-fiber networks under constant stimulation, and has derived de- 
tailed results for several specific types of networks. The following ob- 
servations will treat this same problem, and will presuppose for their 
understanding a reading of the papers mentioned. Our purpose in 
particular will be to develop a somewhat different approach to the 
question; we shall show how this alternative procedure enables us 
both to calculate explicitly all the patterns of steady-state activity 
which are consistent with the applied stimulation in the most general 
possible network—albeit by a rather laborious procedure—and to lay 
the groundwork for a theory of transient excitation. 

We shall make an assumption with regard to the activity of a 
fiber as a function of stimulation which is slightly more general than 
the simple linearity of I, II, III, and IV: we shall suppose that, as the 
stimulation increases, the fiber remains inactive until the threshold is 
reached, then for some range it is a linear function of the excess of 
the stimulation over the threshold; and finally when the stimulation 
reaches a certain upper limit, the resulting excitation has a constant 
value for this and all greater stimulations. This hypothesis seems in 
better agreement with the customary suppositions than that of simple 
linear increase above the threshold, while it occasions us no further 
theoretical difficulties. In particular, the reader will perceive the per- 
durable validity of lemma 1 of Pitts 1942 (hereinafter referred to as 
V) which, indeed, requires little more than that the excitation func- 
tion should be continuous and monotone. 
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In a given network N , which it is convenient to suppose is nota 
simple circuit, we may usefully classify the synapses into three groups. 
The first-order synapses will be those which join links in a simple 
chain, having just one afferent and one efferent fiber ; the second-order 
synapses are simple branch-points, with one afferent and several ef- 
ferent fibers; and the third-order synapses, representing convergence- 
points, possess several afferent and one or more efferent fibers. We 
may simplify N , without changing its activity properties essentially, 
by removing all the second-order synapses; this may be accomplished 
in the following way. If s; be such a synapse, say with n efferent 
fibers d,, do, --- dy, follow the afferent chain to s; contrary to its 
sense until the next preceding second or third-order synapse, say §; , 
is reached. Construct 2 simple chains of fibers such that, for any p, 
the p-th fiber in any of them has the same threshold and activity para- 
meter as the p-th fiber on the chain between s; and s;. Strike out s; 
and the chain from s; to s; , and connect instead the afferent end of all 
the constructed chains to s;, and each one of the n efferent ends to 
one of the 7 fibers which formerly began at s;. Repeat this process 
until no second-order synapses remain. The result will be a network 
N’ which consists of say P third-order synapses, interconnected by 
some Q simple chains. These synapses may be enumerated as s,, 
S.,+++, Sp, and the chains as ¢,, C2,-+-, Cg. Henceforth, we shall use 
the term ‘synapse’ to mean ‘third-order synapse’, when not otherwise 
specified. 


The topological structure of N’, regarded as a complex of chains 
and synapses, may be represented by two functions of integers in the 
following way. ¢(7, 7) = ¢i; is defined to be unity if the chain ¢; 
leads to the synapse s;; otherwise it is zero. y(7, 7) = yi; is unity if 
the chain c; leads from the synapse s; , and otherwise zero. The read- 
er will have remarked that if we consider the chains c¢; as 1-com- 
plexes and the synapses s; as 0-cells in a linear graph, the orientation 
matrix = of N is ||$i;|| — ||yi;||; while conversely, given a 2, we may 
obtain ||¢:;|| and hence ¢ by replacing all negative entries with zero; 
and ||yi;|! by nullifying positive entries and multiplying by —1. 

Various metrical characteristics of the chains ¢; will be impor- 
tant. Among these are the number of fibers in c;, to be denoted by 
‘y;’; the product of the activity parameters of fibers of ¢; , to be desig- 
nated by ‘A,’; and the quantities yu; , defined thus: let « be the exter- 
nal stimulation at the s-th synapse of ¢; (counting its origin as the 
first), less the threshold of the succeeding fibers, and p,, p2, +++ , Dy 
the activity parameters of the fibers of c; , in order: p; is then bs 
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The w; correspond to the quantities o; used for the simple circuit in 
I, Hl, III, and IV. Next, we shall require the quantities y; , which will 
be employed in a slightly different usage than in I, II, III, and IV: y; 
will be the total activity at the synapse s; , diminished by the value of 
o; for this synapse. Further, given a chain c; leading from s; to s;, 
we shall define A; as the greatest lower bound of the values of y; 
which set c¢; into complete activity with the given pattern of external 
stimulation at the synapses within ¢; , and Y; as the least upper bound 
of such values. With respect to each chain c¢;, moreover, and the 
synapse s; whence it leads we may specify a pair of quantities, a; , B:, 
to be called the multipliers belonging to that chain. We shall set 
a; =1if y; 2 Aj; otherwise o, = 0. If y; > Yi, 6; = 0; otherwise 
B; =1. Clearly, its multipliers characterize the state of c; complete- 
ly: in particular, c; is in complete activity if and only if a; 6; =1. 

We may observe first that the state of a network NX over a period 
of time is determined completely by the value of the y; at its third- 
order synapses during that interval; this follows, by a remark made 
in I, from the circumstance that all other synapses are connected to 
a third-order synapse by a simple chain. If, then, we regard the values 
of these y; as components of a point-vector y in a suitable P-space, 
the motion of that point is precisely correlated to changes in the state 
of the network. This representation will be found decidedly useful; 
we may extend it as follows. Along any coordinate axis y; corre- 
sponding to the excitation at s; , we erect perpendicular hyper-planes 
corresponding to the values A;, An, Ar, -*: , Am3 Vz» Yu y*e+) Ym OF 
Yi, Where C;, Cx, ++: , Cm are the chains leading from s;. The totality 
of such planes will divide the whole P-space into boxes or compart- 
ments, each of which will be called a region. Each region I’, corre- 
sponds to a single activity pattern for the chains of N; we shall see 
that the path of the point y, as it moves concomitantly with variation 
in the state of N, has no corners within a region, and only, if at all, 
on the boundaries. 

Thus to every region I’, there is assigned a unique distribution 
of the values zero and one for the multipliers of chains of N , which 
holds for N whenever ¥ is in Ip); and conversely, for every set of 
values for the multipliers which makes no pair of them concomitantly 
zero, we may correlate a unique region. Given I’,, we shall specify 
its multiplier-distribution by the pair of sets 2p,, 2p,, where ap, con- 
tains all multipliers which vanish in Ip, and mp, those which are unity 


there. 
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We may call a point «x of P-space an equilibrium point of NX if 
there exist a steady state of N determined by y = «x and consistent 
with the given external stimulation. A region is stable if it contain 
an equilibrium point; otherwise unstable. It will be seen that the 
steady-state problem treated for various cases in I, II, II and IV is 
essentially a problem in statics, i.e., that of determining all the stable 
regions for a given network having a given applied stimulus pattern. 
This we shall call the static problem; and we shall now proceed to 
find an explicit solution thereof for the general network N. 


II. The Static Network Problem. 

Consider the stimulation at a synapse s; in the steady state. It 
will, in general, be the sum of several types of contribution. First, of 
course, we have the external stimulation o; applied at s;. Second, 
there is the contribution of the chains leading to s;: this may be com- 
puted as follows. Certain of these chains, say ¢;,, Ci,,-+- ; Ci,, Will 
be in complete activity and will consequently deliver an excitation 
Ai, Yi.) + mi,, Where 8,;,; is the synapse where c:, commences. An- 
other group of chains, say ¢;,,, , °°: , Ci,, which are not in complete 
activity ; will deliver generally an excitation Aj, Ai, + ui,; and the re- 
mainder, Ci,,, 5 °° » Ci,, Will contribute an amount A;, Yi,+ i,. 
Adding these, we obtain 


L H K L 
Y=DM, + VA Yu 1g Ai Ai, se Ai Yi. (1) 
G=1 tEal j=H+1 j=K+1 
If we recall some of the quantities defined earlier, we shall be 
able to write (1) more compendiously as 
Yi=Z loi ms + Dvn $76 0 Bs Aj Ye + A; oi [Aj (1 — ay) 
J 
© Xgl = 38s) gs 


where the summations are now taken over all chains and all synapses. 
The apparent extra terms which are in (2) but not in (1) vanish 
through a zero value for some of the multipliers, the ¢:;, or the pj. 

An equation of the kind (2) may be written for every synapse: 
if we put 


M = ||> $4: pix &; Bj Ajl|, 
J 


(2) 


B=|[S¢nls + Apia =o) Ye ea ieee 
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Wwe may write the system in matrix form as 
C—Myj—=R. (4) 


We shall suppose for the moment that the matrix J — M is non-singu- 
lar, and return presently to cases where this does not hold. Under this 
hypothesis, the system (4) may be solved formally for the y; , treat- 
ing the unknown a’s and f’s as indeterminates, by Cramér’s rule: this 
yields a set of expressions for the y; which are the quotients of poly- 
nomials in the multipliers with known coefficients by the determinant 
ie Mi: 
| 
DA; —- M)i;, 5 

where [I — M]j; is the co-factor of the element 7,7 in J — M. 

Now it is clear that any given distribution 1p, , ap, of the multi- 
pliers will determine a stable region if and only if the values for the 
yi resulting by substitution of the given a’s and #’s in (5) satisfy the 
inequalities involved in the definitions of the multipliers. In particu- 
lar, if c; emanate from s;, then the inequalities 

1 > 
i= DRL - Miz = A (6) 
y | fra M| j j 1 
must hold respectively as a ¢ 7%, OF % € ap, , With a similar set for 
the f’s. These requirements may be expressed conjointly by 


(207 =) 4, > (2) 1) Ai 


Um 


a 
(1-2 fy > (1-26) % a) 
for all such j and 1, which is 
> pis (2 a4 — 1) Ry |I a M | ix 
ee ems BESS uae fs) (24, —1) yi; A; > 0 
[2 — M| j 
(8) 


> yi; (1 - 2 B;) Ry | a M\ ix 
ae SD CL = 2:61) Pipl >, 0: 
|I-M | j 
The requirements (8) thus form necessary and sufficient condi- 
tions that a set of a’s and f’s should determine an equilibrium point. 
They may be replaced by a set which is multilinear in the multipliers, 
by making successively the assumptions that |J — M| is positive and 
that it is negative, multiplying by this determinant in (8), and adding 
the corresponding supposition to the set of multilinear inequalities 
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so generated. Higher powers of the o’s and f’s need not occur in 
these, since they may always be reduced to linear occurrences by re- 
membering 0:2 = a; , 6:2 = fi. The resulting conditions will then have 
the form 


s+ P™ a a -ar,¢*+ Bx, Bx, ++ Bx, + Hm > 0. (9) 


LEB Ty eck 
1j,Ky prea ; 


Since the number of multipliers altogether is finite, and each of them 
assumes only the two possible values 0 and 1, the set (9) may always 
be solved in any particular case by substituting in (9) all the possible 
value distributions, and rejecting all which do not fit. There are 3° 
such distributions altogether, Q@ being the number of chains, so that 
the labor might be supposed very great; but in practice the vast ma- 
jority of possible distributions may be found to satisfy or fail to satis- 
fy (9) by inspection. In any case, the problem is one primarily of 
computation; the explicit solutions are not required to deduce the 
general properties of neuron networks, and consequently, while un- 
doubtedly desirable, an expeditious numerical method is not to be re- 
garded as an essential part of the problem; the more so in that discus- 
sions such as our own have as their principal object the enunciation 
of general propositions as to types of possible activity, rather than 
practical application to specific networks. We may remark, inciden- 
tally , that of the distributions 1p, , a, which satisfy (11) , those which 
require both members of some pair a; , 6; to vanish may be discarded 
immediately, as corresponding to no region I’,. This is the source of 
the number 3° given above for the number of possible distributions. 


We may return now to the case where |J — M| vanishes, either 
identically or for certain particular multiplier-distributions. Given 
such a distribution 2p, , t,, we may distinguish two cases. In the first, 
the matrix produced by adding FR as an additional column to J — M, 
(into which we have substituted the values assigned by 7p,, mp, to the 
multipliers) has a greater rank than J — M itself. In this case the 


equations (4) are inconsistent, no vector x of I, satisfies them, and 
there is accordingly no equilibrium point in I’,. The second case, 
where the augmented matrix has the same rank as J — M itself, is 
more interesting. If q be the nullity of I — M, q of the variables y; 
may be chosen arbitrarily, and the equations (4) will then suffice to 
determine values for the remaining coordinates such that the constel- 
lation fixes an equilibrium-point. In this case, therefore, instead of 
having a single equilibrium point in the region, we have a q-dimen- 
sional locus of them, so that even knowledge of both the activity pat- 
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tern and the applied stimulation does not enable us to find the steady- 
state excitation in a unique way. We shall see later, in considering 
the dynamical network-problem, that under these circumstances the 
asymptotic behavior of the system is determined essentially by the 
initial position of the network-vector. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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SOME PROBLEMS IN MATHEMATICAL BIOPHYSICS OF 
VISUAL PERCEPTION AND AESTHETICS 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


The paper outlines a further extension of a previous theory. It dis- 
cusses perception of straight lines, perception of length, perception of in- 
tersections, disconnected contours; effects of associations; breaking up of 
a pattern into parts; some suggestions for experimental procedure. 


In previous publications (Rashevsky 1938, 1940) we have out- 
lined a biophysical theory of visual perception, which among other 
things, led to a physicomathematical theory of aesthetic values of sim- 
ple geometric pattern. The comparison of theoretically calculated 
aesthetic values for 10 simple polygons with experimentally deter- 
mined rank order values for the same patterns gave a rather good 
agreement. Our next step therefore should be to improve and gen- 
eralize the theory, so as to be able to apply it to more complex cases. 
The purpose of the present paper is to outline the way for such a gen- 
eralization. To avoid repetitions we presuppose familiarity of the 
reader with the above two references. 


1. Perception of straight lines. In loc. cit. we assumed as a very 
rough approximation, that the distinct groups of neurons, correspond- 
ing to different vertical components of the eye-movements have all the 
same excitation (Rashevsky, 1938, p. 285; 1940, p. 188). That such a 
simple situation cannot hold exactly, follows from the general discus- 
sion of neural mechanism involved. To each intensity stimulus there 
corresponds a group of neurons lying in different threshold intervals 
h,, h2. In a series of important papers, A. S. Householder (1939, 
1940 a, b) has applied this concept to the theory of intensity dis- 
crimination, which leads to a rather good representation of experi- 
mental data. We may apply his results directly to our problem. 

Using the same notations as Householder (1939; also Rashevsky, 
1940, chap. x), we have for the determination of h, and h, the equa- 
tion: 

AS? ®(Sih) dh 


@(S.h:) =# (Sihe) = (1) 


Met 0 eee OF 
ing 
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The total intensity of excitation of all neurons with thresholds be- 
tween h, and h, is nothing else but the quantity I(S) of Householder, 
and is given by 
AsPa(s, h) dh 
ee 72 

EE) 1 +1(h.— hy) 2) 
If v is the hypothetical constant velocity of scanning (loc. cit.), then 
for the excitation of the V center, the quantity S is given by 


S=S,— av cos &, (3) 


where a is a coefficient of proportionality. For a given ®, which it- 
self is determined by the distribution function f(h) of the thresholds 
of the central neurons, the total intensity of excitation of a V center 
is for every value of the angle 9 given by equations (2) and (3). 

For an H center, S is given by 


Sa Sf bi sinde (4) 


where b is another coefficient of proportionality. Equations (2) and 
(4) determine now the intensity of excitation of the H center for any 
value of 6. 

Assuming different values for the coefficients a and b in equations 
(3) and (4), we ascribe different weights to the vertical and hori- 
zontal components. In this way we should be able to take into account 
the possible preference for the vertical direction, noticed in loc. cit. 
The values of a and b will enter as parameters into the final expres- 
sion for the aesthetic measure of a pattern any may be determined by 
comparison of the theory with experiment. The final aim of the 
theory should be to calculate a and 6 theoretically from assumptions 
about the peripheral neural mechanism. 

As regards the shape of the function 6(S , h), we may try, with 
Householder (1939), the simplest assumption, putting: 


6(S,h) =h(S —h). (5) 


In view of the good agreement of this assumption with data on in- 
tensity discrimination, it will probably suffice as a first approxima- 
tion. Eventually more complex assumptions about 6(S,h) will have 
to be investigated. 

2. Perception of lengths. A similar procedure may be used for 
the calculation of the intensity of excitation of the L center for any 
given length of a straight line. For a very rough approximation we 
may set here S proportional to the length of the straight line per- 
ceived. A recent study of Householder (1940 b) shows, however, the 
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necessity of a somewhat more complex mechanism, and his results 
may be taken over almost directly. 


3. Perception of angles. If the fibers leading to the V or H cen- 
ters send off collaterals to other centers, A*, these collaterals being 
mixed fibers, characterized by either (Rashevsky, 1940) 


A B 
7s ae oe oy Pe af a 5) (6) 
or by 
Amer 


At Be ab: ee (7) 


b ? 
then any sufficiently strong sudden increase of the vertical stimulus 
(8) or of the horizontal stimulus (4) will result in a brief excitation 
of A*, through the collaterals satisfying relations (6). Similarly any 
sufficiently strong sudden decrease of those stimuli will result in an 
excitation of A* through fibers satisfying relations (7). For sake of 
definiteness let us consider the case of a sudden increase of the ver- 
tical component. The other cases are treated in a similar way. 


FIGURE 1 


The stimulus consists now in the sudden increase of the vertical 

component, and is equal (Figure 1) to 
S,, = a,av[cos(6 + ¢ — 2) — cos 6] 
8 
= — a,av[cos(@ + ¢) + cos 6], &8) 

where a, is again a coefficient of proportionality, different from a. 
Thus S,,. is a function both of the magnitude of the angle ¢ and of its 
position, as determined by 0. 

The asymptotic values of « and j at the efferent end of the col- 
lateral fiber of type (6), are, for constant 6 , of the form: 

AE oy BE 


Eo = az Jie baa’ (9) 
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where E is a linear function of S,. The change S,, results in an in- 
crease of FE by an amount AE, AE being proportional to S,,. Hence, 
denoting by a» the coefficient of proportionality, which can readily be 
expressed in terms of the different constants of all the fibers involved 
(Rashevsky, 1938 a, chap. xxii), we have 


AH = 08 Soa (10) 


Because of (6), we have é < Jo. 
The sudden increase AF at the time t = 0 results in a variation 
of « and 7 according to 
AE AAE 
ee Oe 
a a 


oi BE AgBAE ene 
Ls a Sieh eee Gel 


If AE is sufficiently large, then « exceeds 7 at a time ¢, , and then again 
drops below it at the time ¢,. Those two moments are given by the 
roots of the transcendental equation, obtained by putting « = 7 and 
using equations (11). The threshold value of AZ , which must be ex- 
ceeded in order that « shall temporarily exceed 7 is calculated in the 
same manner as the rheobase in the two-factor theory of peripheral 
excitation (Rashevsky, 1938 a). 

A solution of the transcendental equation which determines t, 
and t, cannot be given in a closed form. An immediate practical prob- 
lem therefore is to find a sufficiently good approximation for it. Such 
an approximation gives us ¢, and ¢, in terms of 4H, and hence, be- 
cause of equation (10), also in terms of S,,. Thus 


t, = t (Sp.) ; ty = t2(S,.). (12) 


In loc. cit. we considered, that the temporary excess of « — j at the 
efferent end of the collateral acts as a stimulus upon the afferent end 
of another fiber, whose efferent end is connected to a similar mecha- 
nism as assumed in the theory of intensity discrimination. This will 
result in an excitation of a group of neurons in an “A — center”, such 
that to each positive value of e — 7 at A* (Figure 2) there will corre- 


(v) 


—~ 


\, 
3, 


FIGURE 2 
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spond a distinct group of neurons at A. The relation between « — ] 
at A* and the total intensity of excitation of the corresponding group 
at A is of the same form, as the relation (2), in which ¢ — 7 plays the 
role of S. 

We must however remember, that « — j at A* is not constant, but 
‘varies according to 


=e 9 fade Bait roge0. Bn) 

aaa < iw oil i i ieee itipae ie (13) 
It is positive only between t, and t.. It is effective even within nar- 
rower limits of t¢. If namely h denotes the threshold of fiber IJ (Fig- 
ure 2), then e — 7 is effective only when it exceeds h. This happens 
between two time values t', and t’, > t';, which are the roots of the 
equation in ¢, obtained by equating the right hand side of (13) toh. 
Hence in the expression (2) the intensity 1(S) = I(« — 7) will be 
also a function of time. This function I(t) is obtained by substituting 
in equation (2) « — j for S, and then introducing equation (18) for 
é€— Jj. Since equation (13) depends on AE , and hence on S,, , there- 
fore S,, will also enter as a parameter in the expression for I(t), so 
that we may write /(S,,, t). The total excitation intensity of the A — 
center is given by 


I (Sos) tor = Sp? T (Sos, t) dt. (14) 


By means of equation (8) we may express J(S,.)1.: in terms of 
and @. 
We may also consider an alternate assumption. We may consider 
the value 
Sir(e— jdt, (15) 


which can be obtained from equation (13), as the total value of the 
stimulus S in equation (2). This will give us in general a different 
expression for I(S,.)¢0: in terms of g and @. The technical problem 
involved in both cases is the calculation of t,, t., t,, and t’, in ap- 
proximate closed form. 

The case of a sudden increase of the horizontal component is 
treated in a similar way. For the sudden decrease of either the ver- 
tical or horizontal component we have a very similar problem, except 
that we now havea fiber satisfying relations (7). 

We may consider that all collaterals either of type (6) or of type 
(7), innervated either by the vertical or horizontal component, con- 
verge to the same center A. The total excitation of such a center, 
corresponding to a given angle m, is then obtained by summing the 
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contributions due to S». (increase of vertical component), Sv- (de- 
crease of vertical component); Sy. (increase of horizontal compo- 
nent) and S,- (decrease of horizontal component). 

As an alternative, we may consider however, that each of the 
above four cases corresponds to the excitation of four distinct A cen- 
ters. The perception of an angle would thus always correspond to 
the excitation of four groups of neurons. Analytical expressions for 
the two different assumptions should be derived. 

In loc. cit. we assumed that only the upward scanning of a line 
was effective. In that case, as easily seen, the perception of an angle 
situated as in Figure 3a will consist in complete interruptions of the 
vertical component a cos 6, together with the horizontal component 
b sin 6, (Sy. = — aa cos 6,3 Sy- = — @',b sin 6,), and of the vertical 
component a cos 6, together with the horizontal component 6 sin 42 


FIGURE 3 


(S, = — aa cos 62; Sy- = — @1b sin 6). The total excitation in the 
A center will then depend explicitly only on 6, and 6,, but not on». 
If on the other hand we have the case of Figure 3b, the perception of 
the upper end of the line corresponds only to S,. = — 4@,@ cos 64; 
Su- = a',b cos 6, , which results in a different excitation in the A-cen- 
ter than before. 

For the lower end on the other hand we have S,. = a,a cos 6;, 
Su+ = @1b sin 6,, which gives a still different excitation in the A 
center. 

4. Perception of curves: As outlined in loc cit. the perception 
of a curve involves essentially the same neural mechanism as the per- 
ception of an angle. It is due also to changes in the vertical and hori- 
zontal components, with the difference however, that those changes 
are not sudden but gradual. 

Consider the curve of Figure 4. According to the adopted mech- 
anism (loc. cit.), its perception involves two separate upward scan- 
nings: ABC and ADC. Scanning of AB results in an increase of the 
vertical component from zero to av. This may cause an excitation at 
the efferent end of a collateral of type (6), the excitation being de- 
termined by the shape of the curve. If the equation of AB is given 
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FIGURE 4 


in cartesian coordinates as 
Y=); (16) 
then the distance along the curve of any point A’ from A is equal to 
Pee 
= —) 
F ( 7) x“, (17) 


where x, and «,, are the abscissas of the points A and A’ correspond- 
ingly. Assuming as before, that the curve is scanned with constant 
velocity v , and putting t = 0 for the moment of fixation of the point 


A , we have: 
s=vt. (18) 


Equation (16), (17) and (18) give us x and y for any point A’ of 
the curve as functions of t: 

x=x(t); y=y(e). (13) 
The vertical component stimulus S, at any given point A’ of the curve 
is given by 


df/dx 
Ses 0 COs 0 — a —— (20) 
; V1+ (df/dx)? 
Because of (19), we can express S, as a function of ¢: 
Se eb) (21) 


We thus have a time variable S, , and this will result in a variable ex- 
citation in the V center, J[S,(t)]. The total excitation in the V cen- 
ter, corresponding to the scanning of AB, is given by 


Liot = bee I(S,) dt, (22) 
where #; is the time corresponding to the point B, and is given by 
ta=AB/?. (23) 


Besides the excitation of the V center, a variable S, of the type con- 
sidered will result also in an excitation of the A* center through fibers 
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of type (6). To obtain now the amount of excitation at A*, we re- 
member that at A* 


de dj 
— —=BE — be, (24) 
preter dt y 
where E-is a linear function of S,, and is therefore a known function 
of t. Integrating equation (24) and setting « — 7 = 0, we shall de- 
termine the moments ¢’, and t’, between which « — j > A, and then 
proceed as in section 3; that is either consider 


rts (@ — J) dt 


as the total stimulus, to be used in equation (2) for S, or perform 
the integration over I[(S). 

A general theory of this kind may be of relatively little use, just 
because of its generality. There are, however, definite specific yaue= 
lems that are of practical importance. 

First we may consider AB as an arc of a circle. The differential 
equations (24) are in this case easily integrated, although the deter- 
mination of t’, and ¢’, again requires approximate solutions of trans- 
cendental equations. 

Next we may either study some simplest analytical curves, like 
parabola, cubic, etc. or try to derive some general theorems. One 
might investigate for instance whether a curve whose curvature de- 
creases from A to B gives in general a stronger or a weaker excita- 
tion of the V and A centers, than a curve whose curvature increases 
from AtoB. . 

A scanning of AB results also in a decrease of the horizontal com- 
ponent from bv to zero. In a similar way we may calculate the total 
excitation of the H center during the period of scanning, as well as 
the excitation in the A center, produced through collaterals of type 
Oe 

Similarly scanning BC gives a decreasing S, and an increasing 
S,. The same considerations hold for ADC. We must.keep in mind: 
that for differently variable S, and S, we not only have different in- 
tensities of excitation in the V , H and A centers, but that in each case. 
different groups of neurons are excited. If, as we have done in loc. cit., 
we assume that approximately all those distinct groups have the same 
intensity of excitation, we obtain the following picture for the exci- 
tation pattern of a closed curve, such as shown on Figure 4. Denoting 
by S,., Sy, Sx. and S,- stimuli corresponding to sudden changes due 
to perception of angles, let us denote corresponding gradual changes 


by Sv, Sv-, Sus, Su-. Then the scanning of the curve of Figure 4 
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results in 2 stimuli S,,, 2 of S,, 2 of S,, and 2 of Sy. This gives 4 
groups of neurons in V; 4 groups in H; and 8 groups in A. More- 
over we have 2 groups in the L-center, due to the perception of the 2 
lengths ABC and ADC. Altogether there are 18 groups. Herefrom 
we may calculate the total excitation just as we have done in loc. cit. 
As is easily verified, if we change the shape of the curve, so as to in- 
troduce an inflection point (Figure 5), we obtain a different number 


FIGURE 5 


of groups, since scanning from A to A’ gives first an So , then A’A” 
gives a S,_ and A”B gives again a S,.. With the change of the number 
of neuronic groups the total excitation and therefore the aesthetic 
value of the figure will have changed. 

For a circle we must put, because of symmetry (loc. cit.), the 
two S,. groups as identical, that is as one group with double intensity, 
and similarly for the two S,_ groups. 

With these approximations we could calculate in the manner used 
in loc. cit. the aesthetic values of any combination of straight lines 
and curves. However it is very doubtful, whether the assumption of 
equal intensity of excitation of all groups represent even a rough ap- 
proximation in case of complex patterns. It is therefore more advis- 
able to elaborate the theory of perception of curves in the manner out- 
lined in this section, before proceeding to comparisons with experi- 
ments. 

5. Intersections. Disconnected contours. Consider the point O 
in Figure 6. We actually see four angles at this point. However it 
may be questioned whether the corresponding intensities J in the A 
center will be the same, as if the four angles were lying separately. 
Whereas in the perception of the angle E’, there is no choice and no 
possibility of avoiding the angle, in scanning CEA, in the case of O 


FIGURE 6 
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we may either scan COD, or COA or COB or AOB, ete. For each 
method of scanning we perceive only one angle or no angle at all. 
Therefore, in following such a pattern as in Figure 6, the relative 
frequency of the excitation corresponding to such an angle as AOD 
will be smaller than that corresponding to the angle E'. Accordingly 
we shall attribute a lesser weight to the excitation of the center corre- 
sponding to AOD , by multiplying it by a true fraction, determined by 
the above probability considerations. This fraction is determined by 
the configuration at the intersection; it will for instance be less for 
AOD than for ECO’. 

Similarly we perceive in Figure 6 not only the length CD, but 
also CO’, O'O, OD, CO and O'D. One should investigate whether these 
should be equally weighted. The answer to that question may also be 
found by considering in how many possible ways the whole geometric 
pattern may be scanned and how frequently a particular segment of 
a straight line will be scanned for each way. These questions do not 
arise for simple closed contours, such as discussed in loc. cit. 

Another important problem concerns the neurological mechanism 
of the perception of several distinct, geometrically not connected pat- 
terns. Suppose we see two squares simultaneously. Since we argued 
in loc. cit. that parallel lines, as well as identical angles excite iden- 
tical V , H and A centers, we may think that the total excitation from 
two squares is twice the excitation produced by one. This however is 
not necessarily so, and even appears unlikely. What happens is this: 
In scanning a simple closed contour rapidly, the average intensity of 
excitation of a given center is proportional to the relative frequency 
with which this center is excited per unit time. If we first scan one 
polygon then transfer our eye to another one and scan the latter, the 
number of times a given center becomes excited per unit time remains 
unchanged. An additional excitation is however provided by the move- 
ment of the eye from one polygon to another. This transfer movement 
may itself form a definite pattern. For instance we may have five 
circles, arranged so as to form a regular pentagon. Although, no ac- 
tual pentagon is present in the figure, we perceive this pentagon due 
to the corresponding eye movements when transferring our attention 
from one circle to another. Each movement excites a corresponding 
V , H and perhaps L center. We do not have however the same sensa- 
tion, as when the centers of the circles were actually connected by 
straight lines, forming a pentagon. The intensity of excitation of a 
V or H center is different for an eye movement following an actual 
line, from that when the same eye moves over a blank space. In prin- 
ciple this is not surprising. Besides the excitation of the potential eye 
movement, an actual line produces additional excitation by exciting 
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sO many retinal and central neurons. (Rashevsky, 1938, chap. xxvii). 
The excitation of these neurons may either enhance or inhibit the ex- 
citation of the corresponding V or H centers, depending on the type 
of neuronie connections assumed. By studying different theoretically 
possible cases, we may arrive at a theory of perception and aesthetic 
values of more general patterns, composed partly of imaginary lines. 

6. Effects of associations. It is generally admitted, that any pos- 
sible meaning, associated with a geometric pattern, modifies substan- 
tially its aesthetic value. Neurologically this is also to be expected. 
An “abstract” stimulus pattern produces excitation in a series of cen- 
ters C (Figure 7), which correspond to its geometrical structure. But 


| ! 
ee 


FIGURE 7 


if a pattern reminds us of something, or associates with some other 
previous experience, this means that besides the center C another set 
of centers C’ is excited. If the centers C’ are of an excitatory nature, 
this enhances the total excitation produced by the pattern S, and 
therefore the aesthetic value of S. But C’ may be an inhibitory cen- 
ter, so that its excitation results in an inhibition of other centers. In 
this case our associations reduce the aesthetic value of the pattern S. 
It is however possible that even when C’ is an inhibiting center, the 
association results in an increase of the aesthetic value. This may 
happen, if the ‘associative fibers’, connecting C with C’ are very 
strongly excited and excite, through collaterals, numerous other cen- 
ters. The total excitation is composed of that of C, C’ and the asso- 
ciation fibers. These may be so strong, that the net result is an in- 
crease of total excitation, even though C’ is an inhibitory center. 

These considerations hold quite generally for all types of aesthetic 
experience. We may enjoy a picture that represents realistically strong 
human suffering or we enjoy very sad music (Rashevsky, 1938 a, chap. 
pa m9 TA 

Returning again to the aesthetics of visual patterns, we must re- 
member, that a role similar to that of associations, is played by depth 
and perspective effects. Such effects may also be considered as asso- 
ciational in their nature. Due to the particular arrangement of the 
elements of the pattern, we see them as arranged in three dimensions, 
because they produce associations with centers involved in the per- 
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ception of actual three dimensional figures. 
7” Some general expressions. Consider n centers, of which n, are 
excited with an intensity e,, 2. with an intensity e,, etc. We have 


TuW4=N. (25) 


Let each center inhibit every other one with an intensity j; propor- 


tional to its e; , thus 
Vibe. (26) 


This will be the case, when we neglect the threshold of the inhibitory 
fibers (Rashevsky, 1938 a, chap. xxii). The more general case of 
finite thresholds is treated in a similar way. 

Each center in the group 1 receives an amount of inhibition: 
b, (1, — 1) e; + benzez + -+- =D njbie; — b,e, and more generally , each 


center of the group k receives an amount of inhibition equal to 


D> nidbie; — Oxy. 
i 


Since there are m, centers in the k-th group, the total inhibition in that 
group is 2% (> 2:bie; — b,e,). The total inhibition is obtained by sum- 
ming this with respect to k , and because of (25) is equal to 


j= (i 1) ben, (27) 
The net excitation is equal to 
ES ne; =] (a1) > bes. (28) 
If all b; have the same value b, equation (28) reduces to 
H=[1— (n=1) b] Sine;. (29) 


This is the case, considered in loc. cit. 

When 7, which measures the complexity of the pattern, exceeds 
(1 + b)/b, then E becomes negative. This however does not mean 
that all centers are inhibited, for EF gives the net excitation. Some 
centers remain excited, some are inhibited, even when E < 0. (Ra- 
shevsky, 1938a, chap. xxii). On the other hand, some centers may be 
already inhibited even when n < (1+ 5b)/b, so that the net value of 
E is still positive (Rashevsky, 1938 a, chap. xxii). 

The excited centers are those, for which e > Snjbje; — bree . 
Since all other centers are inhibited, we arrive at first glance at the 
absurd conclusion, that some elements of the design, when the latter 
is too complex, will not be perceived at all. This paradox is resolved 
however in the following manner: Our attention is concentrated al- 
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ternately first on those elements of the design, which correspond to 
the excited centers, then on those which correspond to the inhibited 
ones. But in this latter case we have to take in the expression for E 
not the total x, but only the value n* corresponding to the inhibited 
elements. This makes n smaller, and reduces the complexity, so that 
all centers now may be excited. In other words, we now break up the 
very complex pattern in two (or more) partial patterns of lesser com- 
plexity, and we either perceive the total pattern as consisting of two 
or several distinct parts , or as consisting of one outstanding part, 
seen on a general background. 

Such a splitting of a very complex pattern is actually observed 
on some of the patterns, used by L. Wiley (1940) in his work on eval- 
uation of aesthetic values of abstract designs. For instance in his de- 
sign No. 9 we definitely see the two block-like figures on a general 
background of criss-cross lines. Until we have an adequate theory, 
developed along lines suggested in preceding sections, we cannot pre- 
dict exactly what part of the design should stand out. Using the crude 
approximation, adopted in loc. cit. and discussed for curves on p. 185 
of this paper, however, we find that the strongest excited centers, 
which would remain uninhibited even if the pattern would be per- 
ceived at once, and not in parts, are those corresponding to the lines 
and angles forming the block-like figure which is seen to stand out on 
the background of the remaning elements. A similar situation holds 
for a few other of the complex designs, used by Wiley. 

8. Some suggestions for experimental procedure. The final goal 
of a biophysical theory of visual aesthetics is the theoretical calcula- 
tions of the aesthetic values of geometric patterns of any complexity 
and its comparison with psychophysically determined rank order val- 
ues. However, this goal cannot be reached before all the preliminary 
work, outlined in this paper has been done. As this work on the par- 
tial problems proceeds, experimental verifications should be done all 
along. 

The simplest possible way for such preliminary experiments 
would be to study first straight lines alone, then angles alone, etc. This 
cannot be done however, because we cannot perceive an angle without 
having at least two straight lines to form it. Neither can we have a 
piece of straight line without ends, the perception of which, as we 
have seen in section 3, involves some A-centers, the centers involved 
varying with the position of the straight line in the plane. A com- 
promise method must therefore be adopted. | 

Having developed a more exact theory of the perception of an- 
gles we may verify it by experimenting on different designs, in which 
lengths and positions of straight lines in the plane do not vary, but 
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angles do. Such is the case with Figures 8a and 8b. For the study of 


{} 


A 8 c 


FIGURE 8 


the effects of the L-center we may compare designs like those of Fig- 
ures 8a and 8c, as well as similar polygons of different sizes. In con- 
sidering similar curvilinear designs, we must remember that an in- 
crease in size means a decrease of curvature, and that this may in 
general affect the intensity in the A-centers involved. Having thus 
definitely established different quantitative aesthetic relations pertain- 
ing to angles and lengths, we can use designs involving different V 
and H centers.. Even if the L and A centers now vary from design to 
design, their effect will be known from the preceding study. 

A comparison of aesthetic values of designs, consisting of poly- 
gons drawn in full lines and similar polygons indicated by dots at 
their vertices or by some simple pattern, arranged at the vertices, will 
contribute to the studies discussed in section 6. 

As has been stated in loc. cit. different parameters entering into 
our equations and determining the theoretical rank order for a given 
series of design, vary from individual to individual. This takes into 
account individual tastes. In actual rank order measurements we de- 
termine some sort of average values. Those average values will de- 
pend on the distribution function of the different parameters. A theo- 
retical study of the effects of different distributions will form an im- 
portant problem. With the aid of such studies important conclusions 
may be drawn from the distribution of individual scores for a given 
design. 

A complete mathematical theory would also indicate to what ex- 
tent the total aesthetic measure of a design varies with any of the 
parameters. If a small variation of one of the parameters may result 
in a large variation of the aesthetic value of the design, we should ex- 
periment with designs that are sufficiently different in their aesthetic 
values. A difference of experimental rank orders 8-5 and 8-4 may be of 
no use, and either an agreement or disagreement of theory and ex- 
periment in such a case may be meaningless. Thus it may be better 
to use a fewer number of designs spread sufficiently apart, rather than 
a very large number, forming an almost continuous series. This will 
also result in a considerable saving of time and labor. 
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